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Nearby supernova remnants and the cosmic-ray spectral hardening 
at high energies 
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ABSTRACT 

Recent measurements of cosmic-ray spectra of several individual nuclear species by the 
CREAM, TRACER, and ATIC experiments indicate a change in the spectral index of the 
power laws at TeV energies. Possible explanations among others include non linear diffusive 
shock acceleration of cosmic-rays, different cosmic -ray propagation properties at higher and 
lower energies in the Galaxy and the presence of nearby sources. In this paper, we show that 
if supernova remnants are the main sources of cosmic rays in our Galaxy, the effect of the 
nearby remnants can be responsible for the observed spectral changes. Using a rigidity depen- 
dent escape of cosmic-rays from the supernova remnants, we explain the apparent observed 
property that the hardening of the helium spectrum occurs at relatively lower energies as com- 
pared to the protons and also that the spectral hardening does not persist beyond ~ (20 — 30) 
TeV energies. 
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1 INTRODUCTION 

Recently, cosmic-ray (CR) measurements by the new-generation 
balloon-borne experiments such as the ATIC (Panov et al. 2007), 
CREAM (Yoon et al. 201 1), and TRACER (Ave et al. 2008) seem 
to indicate that the CR spectrum deviates from a single power law. 
The spectra of all individual elements seem to be harder at TeV 
energies than at lower energies. Such a hardening is not easy to 
explain under the standard models of CR acceleration and their 
propagation in the Galaxy. Under the standard theory, CRs below 
the knee (~ 3 PeV) are considered to be produced by supernova 
remnant (SNR) shock waves by diffusive shock acceleration mech- 
anism (Bell 1978, Blandford & Eichler 1987). Such a mechanism 
naturally predicts a power law spectrum of E ' with the index 
7 = 2 for strong shocks. On the other hand, CR propagation in 
the Galaxy is considered to be of diffusive nature which is due to 
scattering by magnetic field irregularities and the CR self excited 
Alfven and hydromagnetic waves present in the Galaxy. Measure- 
ments of CR secondary-to-primary ratios indicate that the diffusion 
is energy dependent with the diffusion coefficient D(E) oc E a 
with a » (0.3 — 0.7). Under these considerations, the CR spec- 
trum in the Galaxy is expected to follow a single power law with 
index (7 + a) up to the knee, which do not seem to agree quite 
easily with the observed hardening at TeV energies. 

The observed data can be explained if either the source spec- 
trum or the diffusion index flattens at higher energies. Non linear 
diffusive shock acceleration theories where CRs modify the shock 
structure predict concave spectra (flatter at higher energies) at the 
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shocks. But, the total spectrum injected into the interstellar medium 
which is the sum of the instantaneous spectra over the SNR life time 
is very close to a pure power-law (Caprioli et al. 2010). The concave 
signature can be even more diluted when summed over an ensem- 
ble of SNRs (Ptuskin & Zirakashvili 2005). From the propagation 
point of view, there are models which assume a harder or constant 
CR diffusion coefficient at higher energies in the Galaxy (Ave et 
al. 2009). Such models are motivated not only by the apparent flat- 
tening of the observed boron to carbon ratio above ~ 100 GeV 
energies, but also by the observed CR anisotropy which is almost 
independent of energy. Recently, it has also been proposed that dis- 
persion in the spectral indices of CR source spectrum from many 
sources can also be responsible for the observed spectral hardening 
(Yuanetal. 2011). 

Another possible explanation, as also pointed out in Ahn et 
al. 2010, is the presence of nearby sources. Erlykin & Wolfendale 
201 1 suggested that an extra component of CRs with a steep spec- 
trum could be contributing below ~ 200 GeV/n while above that, 
the spectrum is entirely determined by a harder CR background. 
They proposed that the sources of the extra component could be in 
OB associations in the Local Bubble. Recently, Ohira & Ioka 201 1 
proposed that the hardening could be due to decreasing Mach num- 
ber in hot superbubbles with multiple supernovae. In another recent 
work, Vladimirov et al. 201 1 investigated several possible interpre- 
tations (including local source effect) for the observed spectral fea- 
tures at low and high energies using the GALPROP propagation 
code. They also presented the possible effects on other observed 
properties such as CR anisotropy, isotopic ratios and the Galactic 
diffuse 7-ray emissions. 

In our present study, we investigate whether the spectral hard- 
ening observed at TeV energies could be an effect of the nearby 
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SNRs. Although there has not been any direct detection of CRs 
from any sources, SNRs remain the most favorable candidates both 
theoretically as well as observationally. At least the presence of 
high energy particles up to few TeVs inside SNRs have been con- 
firmed by the detections of non-thermal X-rays (Parizot et al. 2006) 
and TeV 7-rays from several SNRs (Aharonian et al. 2006, 2008a). 
Moreover, the detection of TeV electrons by the HESS experiment 
(Aharonian et al. 2008b) indicates the presence of one or more CR 
sources within a distance of ~ 1 kpc from us. If these sources pro- 
duce both, electrons as well as nuclei, we expect to see some effects 
on the spectra of CR nuclei observed at the Earth. 



2 MODEL 

The diffusive propagation of CRs in the Galaxy neglecting the ef- 
fects due to nuclear spallation can be described by the following 
equation, 



dN 

V-(DVN)+Q=^ 



(1) 



where N(r, E, t) is the differential number density at a distance 
r from the source at time t, E is the kinetic energy/nucleon and 
Q(r,E,t) is the source term. The diffusion coefficient is taken 
as L>(K) = D (K/5fto) a for 5ft > » , where K denotes the par- 
ticle rigidity which is given by 5ft = AE/Z for charge Z and 
mass number A. For our study, we consider two sets of values 
for (Do , Sfto , a) : one based on purely diffusion model (hereafter 
Model A) and the other based on models including CR reaccelera- 
tion due to interstellar turbulence (hereafter Model B). We choose 
(D ,% h a) = (2.9,3,0.6) for Model A (Thoudam 2008) and 
(5.8, 4, 0.33) for Model B (Strong et al. 2010), where D is in 
units of 10 2S cm 2 s" 1 and K is in GV. 

Under diffusive shock acceleration theory, CRs are confined 
within the remnant due to the magnetic turbulence generated by 
the CRs themselves. They can escape when their upstream diffu- 
sion length defined as luff = D s (E) /u s is greater than the escape 
length from the shock front which is usually taken as l eac w 0. 1R S , 
where u s and R s denote the shock velocity and the shock radius re- 
spectively. In the Bohm diffusion limit, the upstream diffusion co- 
efficient scales linearly with energy as D s (E) cc E which implies 
that higher energy particles can escape the remnant at early times 
followed later by the lower energy ones. But, the exact energy de- 
pendence of D s is still not well understood and depends on some 
poorly known quantities which include the spectral distribution of 
the CR self-excited turbulence waves, the level of magnetic field 
amplification by the CRs and the dynamical reaction of CRs on the 
shock structure. Therefore, we follow a simple but reasonable pa- 
rameterization for the CR escape time similar to that adopted by 
Gabici et al. 2009 as given below, 



^esc(3^) — ised 



-1/a 
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where t se d denotes the start of the Sedov phase, K maa; denotes 
the maximum CR rigidity and a is a positive constant. We assume 
that the maximum CR energy accelerate by an SNR scales with the 
charge number Z as ZU max , where U max denotes the maximum 
kinetic energy of the protons which is taken as 1 PeV for our study 
(Berezhko 1996). This scaling gives K ma ^ = 1 PV. In units of 
energy/nucleon, the maximum energy for helium is E max = 0.5 
PeV/n. 

Eq. (2) assumes that the highest energy CRs of all the species 



start escaping at the onset of the Sedov phase itself. Writing Eq. 
(2) in terms of total kinetic energy, it is easy to check that for the 
same kinetic energy, the escape time of CRs scales with the charge 
number as Z 1/a , i.e, higher charged particles escape at relatively 
later stages of the SNR evolution. Thus, our escape model takes into 
account the general understanding of diffusive shock acceleration 
theory that higher charged particles can be confined for relatively 
longer duration within the remnant. In terms of energy/nucleon, we 
can write Eq. (2) as 
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Eq. (3) shows that for the same energy/nucleon, all nuclei with 
charge Z > 1 escape earlier than the protons by a factor of 
(A/Zy 1/a . We further assume that no particles remain confined 
after the shock completely dies out which we assume to occur when 
the SNR age 10 5 yr. Taking this into account, the CR escape time 
for our study is taken as T esc (E) — min \t esc (E), 10 5 yr] . For 
detailed studies on particle escape from SNRs, see e.g., Ptuskin & 
Zirakashvili 2005, Caprioli et al. 2009 and Ohira et al. 2010. 

The corresponding escape radius of CRs is calculated using 
the age-radius Sedov relation for SNRs as given below, 



Resc(E) = 2.5ti t sea 



-0.6 



(4) 



where uo represents the initial shock velocity, i.e the velocity at 
The source term in Eq. (1) is taken as, 



Q(r,E,t) 



9(E) 



S(t - T esc )5{r - R e3C ) 



(5) 



where A esc — 4ttR 2 sc denote the surface area of the SNR at the 
time when CRs of energy E escape the remnant. It should be noted 
that our consideration of the rigidity dependent escape time and the 
finite source size are different from the commonly adopted burst- 
like point source approximation where CRs of all rigidities are as- 
sumed to escape at the same time from a point source. For CR 
study near the sources, the point source approximation can break 
down and it looks more realistic to take their sizes into account 
(Thoudam & Horandel 2011). Recently, such importance has also 
been highlighted in Ohira et al. 201 1 in the study of 7-ray emission 
from SNRs interacting with molecular clouds. 

The source spectrum in Eq. (5) is taken as q{E) = Aq(U) 
with q(U) given by, 



q(U) = k{U 2 + 2Um)- ( '< +1)/2 {U + rn) 



(6) 



where U = AE represents the particle total kinetic energy, m is 
the rest mass energy and k is the normalization constant which is 
related to the CR injection efficiency. 

Solving Eq. (1), the spectrum at a distance r s from the SNR 
follows, 



N(r 3 ,E,t) 



q(E) R e 



1 s-f^esc 



=exp 



T e3C j 
x sinh 



R 2 + r 2 ) 



AD(t - Tesc 
T s Resc 



K 2D(t-T esc ) 

For high energy particles for which the diffusion radius defined as 
rdiff = \/D(t — Tesc) is much larger than (r s ,R esc ), Eq. (7) 
follows a power-law of the form N(E) oc E~( r+ % a ) . 

Eq. (7) can be used to calculate the CR spectra from the nearby 
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SNRs. We choose proton and helium for our study and consider 
only those SNRs with distances < 1 kpc from the Earth and ages < 
2x10 yr. From the available literature, we found 10 SNRs listed as 
follows with their distances (kpc) and ages (yr) respectively given 
in parentheses: Cygnus Loop (0.54, 10 4 ), HB21 (0.8, 1.9 x 10 4 ), 
HB9 (0.8, 6.6 x 10 3 ), S147 (0.8, 4.6 x 10 3 ), Vela (0.3, 1.1 x 10 4 ), 
G299.2-2.9 (0.5,5 x 10 3 ), SN185 (0.95, 1.8 x 10 3 ), Monogem 
(0.3,1.1 x 10 s ), G114.3+0.3 (0.7,4.1 x 10 4 ) and Vela Junior 
(0.75,3.5 x 10 3 ). 

In addition to the contributions from the nearby SNRs, we as- 
sume that there exists a steady CR background in the Galaxy which 
dominates the overall CR spectrum. For the CRs observed at the 
Earth, we assume that this background component consists of con- 
tributions from distant SNRs plus any other possible sources in the 
Galaxy. For our study, we obtain the background by fitting the ob- 
served CR spectrum between (20 — 200) GeV/n. This is the en- 
ergy region where the contamination from the nearby sources is 
expected to be less and at the same time, not much affected by the 
Solar modulation. In fact, it has been shown in Thoudam 2008 that 
the presence of nearby sources can produce stronger density fluctu- 
ations at higher energies than at lower energies because of the en- 
ergy dependent nature of CR diffusion. Therefore, we believe that 
it is reasonable to assume that the low energy CRs that we observe 
at the Earth are not much affected by the presence of nearby SNRs 
and they largely represent the averaged background spectrum in the 
Galaxy. We will show in the following that this is indeed the most 
likely case. 



3 RESULTS 

From the fit, the spectral indices of the background CRs are found 
to be 2.75 ± 0.01 for the protons and 2.68 ± 0.02 for the helium. 
The reason for the flatter helium spectrum is not properly under- 
stood. Recently, Blasi & Amato 2011 showed that the flatter he- 
lium spectrum with respect to the protons above 1 TeV could be 
due to spallation effects. Later, Vladimirov et al. 201 1 showed that 
such effects can lead to boron-to-carbon ratios and the anti-proton 
fluxes which are inconsistent with the observed data. Another pos- 
sibility for the different spectral indices could be that the intrinsic 
source spectra itself are different. It could be due to different ac- 
celeration sites of protons and helium (Biermann et al. 2010) or 
inhomogeneous abundance of elements in superbubbles (Ohira & 
Ioka 201 1). For our present study, we assume that CRs are injected 
into the Galaxy with different source indices. The index 7 for an in- 
dividual species is chosen such that (7 + a) is equal to the spectral 
index of the background obtain from the fit. 

Before illustrating our results, we briefly discuss the choice 
of other model parameters involved in our calculations. Typically, 
t B ed has values between ~ (100 — 10 3 ) yr depending on the gas 
density of the interstellar medium, mass of the ejecta and the energy 
output of the supernova explosion. For our study, we take t se d = 
500 yr. We assume the initial shock velocity uo to be 10 9 cm/s. 
This gives CR escape times from the SNRs in the range of t eac = 
(500 — 10 J ) yr and the corresponding escape radii as R eS c ~ (5 — 
100) pc. Finally, we treat the escape parameter a and the injection 
efficiency of the protons (helium) hereafter denoted by £ p (h e -s as 
free parameters. For our calculations, we will assume that all the 
parameters mentioned above are same for all the SNRs. 

Because of lack of precise informations on the values of a and 
£p(he) > we perform calculations for several of their randomly cho- 
sen combinations. We choose the escape parameter in the range of 



CREAM 
ATIC 
AMS 
BESS 
PAMELA 

Background: Fitted 
Random spectra 
Averaged spectrum 



Proton (Model A) 




10 s 10* 
E (GeV) 



CREAM 
ATIC 
AMS 
BESS 
PAMELA 

Background: Fitted 
Random spectra 
spectrum 



Proton (Model B) 




E (GeV) 



Figure 1. Proton spectra (xE 2 - 75 ) for Model A (top panel) and Model B 
(bottom panel). The blue dotted line represents the background spectrum. 
The thin black lines represent an example of 30 random spectra calculated 
with proton escape parameters and the injection efficiencies in the range of 
a = (1 — 3) and e p = (5 — 25) x 10 49 ergs respectively. Each spectra 
is the sum of the background and the contribution from the nearby SNRs. 
The blue solid line represents the average of 200 random spectra. See text 
for data references and other details. 



a = (1 — 3). This range approximately covers the a values given 
in some available literatures. Studies based on non-linear diffusive 
shock acceleration theories which takes into account the modifica- 
tion of the shock structure by the CRs give a ~ 0.8 (e.g., Ptuskin & 
Zirakashvili 2005). Blasi & Amato 2011 adopted a ~ 3.2 in their 
study of the effect of random nature of SNRs on the CR spectrum. 
Investigations of 7-ray emissions from molecular clouds interact- 
ing with nearby SNRs adopt values in the range of a = (2.4 — 2.6) 
(Gabici et al. 2009, Ohira et al. 2011). We consider the CR injec- 
tion efficiencies in the range of e p — (5 — 25)% for protons and 
£he = (1 — 5)% for helium, where the values are in units of 10 51 
ergs. The averaged proton to helium injection ratio of 5 which we 
consider here is less than the observed proton to helium flux ratio 
of ~ (20 - 13) in the energy range of ~ (20 - 200) GeV/n (Yoon 
et al. 201 1). But, our wide range of efficiencies for both the species 
well cover the observed flux ratios. It should be understood that the 
observed flux ratios may not necessarily represent the injection effi- 
ciency ratios from the source. Effects during the propagation in the 
Galaxy such as due to spallation (which are different for different 
nuclear species depending on their interaction cross-sections) may 
change the composition ratios produced by the source. In addition, 
propagation of CRs is charged dependent. Those which undergo 
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Figure 2. Helium spectra (x_E 2,7 ) for Model A (top panel) and Model B 
(bottom panel). The calculation assumes injection efficiencies in the range 
of e he = (1 — 5) X 10 49 ergs. All other model parameters, results repre- 
sentation and the data references remain the same as in Fig. 1 . 



faster diffusion will escape more easily from the Galaxy and even- 
tually lead to less flux observed at the Earth. 

Fig. 1 shows our calculated proton spectra (xi? 2,75 ) for Mod- 
els A (top panel) and B (bottom panel). In each panel, the thin 
black lines represent an example of 30 different random spectra 
we have calculated. Each random spectrum corresponds to a set of 
{a, t v {he.) } which is assumed to be the same for all the SNRs. Each 
spectrum is the sum of the background CRs (shown as the blue dot- 
ted line) and the total contribution from all the nearby SNRs we 
have considered. The blue solid line represents the averaged spec- 
trum of a total of 200 such random spectra. The data are taken from 
CREAM (Yoon et al. 201 1), ATIcfl (Panov et al. 2007), AMS (Al- 
caraz et al. 2000, Aguilar et al. 2002), BESS (Haino et al. 2004), 
and PAMELA 1 (Adriani et al. 2011) experiments. One common 
result that we can notice between the two models is that the con- 
tribution of the nearby SNRs show up mostly above ~ (0.5 — 1) 
TeV. However, there are some general differences between the two 
results. The results for Model A not only show larger variations 
between individual spectra but also stronger irregular features and 
spikes. Also, in general Model A produces larger contribution from 
the nearby SNRs as compared to Model B. This is largely due to the 
comparatively harder source spectrum of CRs required in Model 
A. For the reasonable range of injection efficiencies considered in 



1 Data taken from the database compiled by Andrew W. Strong (Strong & 
Moskalenko 2009) 



our study, the results of Model A seem to be in better agreement 
with the data both in terms of the size and the shape of the spec- 
tra. On comparing the averaged spectra (thick blue lines) above 
~ (0.5 — 1) TeV, the result of Model A is comparatively harder up 
to ~ 100 TeV which then becomes steeper at higher energies. This 
spectral behavior of Model A is in good agreement with the recent 
data which also seem to indicate that the spectral hardening for pro- 
tons does not persist beyond ~ (20 — 30) TeV. On the other hand, 
the averaged spectrum in Model B show less hardening above ~ 1 
TeV and it continues without any turn over or steepening up to the 
maximum energy considered here. 

The corresponding results for helium are shown in Fig. 2: 
Model A (top panel) and Model B (bottom panel). Our results for 
helium look similar to those obtained for protons. One general dif- 
ference we notice is the shifting of the helium results towards lower 
energies with respect to the proton results. Though not very signif- 
icant, a similar trend is also present in the observed data. For in- 
stance, the spectral hardening in the helium data occurs at ~ 0.5 
TeV/n whereas for the protons it occurs at ~ 1 TeV. Moreover, the 
spectral turnover at higher energies seems to occur at ~ 10 TeV/n 
for helium while for protons it seems to occur at ~ (20 — 30) TeV. 

In Fig. 3, we present our best fit results: protons (xE 275 , 
top panel) and helium (xE 2 ' 7 , bottom panel). They are obtained 
by choosing (a,e p ,e h e) = (2.2, 9%, 2%) for Model A and 
(2.4, 20%, 3.7%) for Model B. Our model parameters give escape 
times of t esc — (500 — 10 ) yr for protons of energies (1 PeV— 8.6 
GeV) and for helium of (0.5 PeV/n-4.3 GeV/n) in Model A. In 
model B, the corresponding values are (1 PeV— 3 GeV) and (0.5 
PeV/n— 1.5 GeV/n) respectively. The data in Fig. 3 are the same as 
in Figs. 1 & 2 respectively. The blue dotted line represents the back- 
ground CR spectrum. The solid lines correspond to Model A and 
the double dotted lines to Model B in which the thin and the thick 
lines represent the total contributions from the nearby local SNRs 
and the total background plus nearby contributions respectively. In 
Model A, the dominant nearby contributors are the Vela, G299.2- 
2.9 and SN185 remnants. They are shown by the thin dashed lines 
in the figures. Vela dominates in the range of ~ (0.7 — 10) TeV/n 
while above that, the spectrum is mostly dominated by G299.2-2.9 
and SN185. In Model B, Vela dominates over a wide range up to 
~ 300 TeV/n and beyond that, it is dominated by G299.2-2.9 (not 
shown in the figure). 

The steep low energy cut-offs in the individual SNR contribu- 
tions in our model are largely due to the energy dependent escape 
of CRs. CRs below the cut-offs are mostly those which are still 
confined within the SNRs and are not yet released into the inter- 
stellar medium. Our best fit result for Model A show a rise in the 
total spectrum at ~ (0.5 — 1) TeV/n which remain constant up to 
~ (5 — 10) TeV/n. This is due to the effect of the low energy cut-off 
of the Vela remnant. To be specific, the rise in the proton spectrum 
occurs at ~ 1.2 TeV while that of the helium at ~ 0.6 TeV/n. 
This difference is largely due to the effect of early escape of he- 
lium compared to the protons for the same energy/nucleon. There 
is also some effect due to the faster diffusion of helium than the 
protons for the same energy/nucleon. As mentioned in section 3.1, 
it is interesting to see that the recent data also seem to indicate that 
the helium spectrum starts hardening at comparatively lower ener- 
gies than the protons. Our best fit results then show a slow increase 
above ~ (5 — 10) TeV/n which again becomes almost constant 
above ~ (40 — 60) TeV/n. This is due to the combined effect of 
other SNRs mainly G299.2-2.9 and SN185. These spectral features 
are found to be more pronounced for helium. In Model B, these 
features are smeared out and we get a smooth spectrum with a slow 
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Figure 3. Best fit results for proton (xE 2 - 75 , top panel) and helium 
(xE 2 - 7 , bottom panel). The data are the same as given in Figs. 1&2 re- 
spectively. The blue dotted line represents the background spectrum. Solid 
lines correspond to Model A and double dotted lines to Model B of which 
the thin line represents the total contributions from the nearby SNRs and 
the thick line represents the total background plus nearby contributions. The 
thin dashed lines represent the dominant nearby contributors in Model A. 



increase above a few TeVs. This is due to the comparatively slower 
diffusion of CRs in this model and the dominance by a single source 
(Vela) over a wide range of energy spectrum. 

The total CR anisotropy A expected under our model can be 
calculated using the following equation (Thoudam 2007), 



A 



(8) 



where the summation is over the nearby SNRs, f j denotes the di- 
rection of the i th SNR giving an intensity Ii at the Earth, f m de- 
notes the direction of maximum intensity, It represents the total 
observed CR intensity and Si denotes the anisotropy amplitude due 
to a single SNR. Si under the diffusion approximation is given by 
(Mao & Shen 1972), 



Si = 



3D \VNj 
c N t 



(9) 



denotes the CR density from an SNR 
For our best fit proton results, we get 



where Ni (given by Eq. 7) 
with distance r, and age t 4 . 
A ss (1.7 x 10~ 2 - 0.12) and (1 - 4) x 10~ 2 for Model A and 
B respectively in the energy ] 
are larger than the measured ; 
the same energy range. But, 
more closer to the measured ' 



■ range of (1 — 100) TeV. Our estimates 
I anisotropics of ~ (0.5 — 1) X 10 -3 in 
, compared to Model A, Model B looks 
i values (see also Ptuskin et al. 2006). 



4 DISCUSSIONS AND CONCLUSIONS 

We show that for both Models A and B, the nearby SNRs contribute 
mostly above ~ (0.5 — 1) TeV/n and they may account for the ob- 
served spectral hardening at high energies. We show this for a wide 
range of CR injection efficiencies and CR escape parameters from 
the SNRs. Looking into the averaged spectra in Figs. 1&2, we find 
that both the models predict that the hardening of the helium spec- 
trum should occur at lower energies as compared to the protons. 
We also find that the averaged result of Model A seems to explain 
the overall data better than that of Model B. 

However, the wide range of parameters values considered in 
our study allow both the propagation models to successfully ex- 
plain the observed data with a careful choice of model parameters. 
We show this with our best fit results in Fig. 3. But, the high CR 
injection efficiency of e p — 20% required in Model B is around 
a factor of 2 larger than the normally considered value of ~ 10% 
for CR studies in the Galaxy. Moreover, the steep source index of 
7 ~ 2.4 required in this model is also hard to reconcile with the 
results of diffusive shock acceleration theory which predict an in- 
dex of 7 = 2. Model A, on the other hand, looks favorable con- 
sidering its relatively more reasonable values of the source index 
(7 = 2.15) and the proton injection efficiency (e p = 9%) required 
to explain the observed hardening. In addition, Model A also bet- 
ter explain the apparent observed property that the spectral hard- 
ening does not persist above a few TeVs. However, the measured 
anisotropy seem to favor Model B which assumes a weaker energy 
dependence of CR diffusion in the Galaxy. 

Our results look different from the predictions of other mod- 
els. Models based on constant diffusion coefficient at high energies 
or spectral dispersion in the source spectrum are expected to pro- 
duce a high energy spectrum which remains hard up to the maxi- 
mum energy (Ave et al. 2009). But, the data indicates that the spec- 
tral hardening happens only up to ~ (20 — 30) TeV for protons 
and ~ 10 TeV/n for helium which in general agrees well with our 
predictions (especially with Model A). It should be mentioned that 
our results may not be significantly different from others if the CR 
spectrum has a break or a cut-off (normally assumed to be expo- 
nential) at energies < 0.1 PeV. In such a case, the spectrum will 
start rolling over before it starts showing noticeable differences. 
But, note that a cut-off somewhere between ~ (3 — 5) PeV is pre- 
ferred, irrespective of the nature and the origin of the cut-off, in 
order to explain the observed knee in the energy spectrum of CRs 
(Horandel 2003). 

The secondary CR spectrum under our model can be even 
more different from other models. Secondaries are those which 
are considered to be produced by the spallation of the primaries 
only during the propagation in the Galaxy. Their spectrum 7V S in 
the Galaxy is related to their primary spectrum N p as N S (E) oc 
N P (E)/D(E). Thus, for N P (E) oc E~ T , the secondary spec- 
trum follows N S (E) oc _E _ ( r+a ) which is steeper than their pri- 
maries by the diffusion index a. Therefore, once D(E) is fixed, 
the shape of the secondary spectrum depends on the shape of their 
primary spectrum. This means that models which assume the same 
D(E) but different N P (E) will produce different N S (E). Under 
our model, if we neglect the production of secondaries inside the 
SNRs, we can assume that all the secondaries are produced by the 
background CRs. As our background primary spectrum is steeper 
than the spectrum used in other models to explain the spectral hard- 
ening, e.g. Yuan et al. 2011, we expect a steeper secondary spec- 
trum in our case. This difference would be even more significant 
when compared to propagation models which assume a constant 
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CR escape time from the Galaxy at higher energies (Ave et al. 
2009). Under such models, N S (E) oc E~ r at higher energies while 
at lower energies N B (E) oc E~ (r+a) . The differences we just 
mentioned are expected in all kinds of secondary nuclear species 
like boron, sub-Fe, and anti-protons. At present, data on secondary 
spectra are available at most only up to ~ 100 GeV/n. Future high 
energy measurements would be crucial to test our model. 

In addition, the diffuse 7-ray emission of our Galaxy can also 
provide an important check of our model. If the diffuse emission 
is dominated by the 7r°-decay 7-rays, then their intensity would 
largely follow the proton spectrum at high energies. Therefore, 
under our model we expect a diffuse spectrum which is steeper 
than the predictions from other models. In fact, it has already been 
shown in Yuan et al. 2011 that under their model, the 7-ray spec- 
trum would become harder above ~ 50 GeV. Preliminary results 
from the FERMI measurements up to ~ 100 GeV show that the 
spectrum is in good agreement with the conventional model as- 
suming a single power-law CR spectrum above a few GeV (Strong 
201 1). The spectrum do show some excess above the model which 
could well be attributed to unresolved point sources like pulsars. 
Detailed investigation of the diffuse 7-ray spectrum and also future 
measurements at even higher energies would be important to check 
the validity of our model. 

In short, we conclude that the apparent change in the spectral 
index of the CR energy spectra at TeV energies could be a local 
effect due to nearby SNRs. A detailed investigation of both the 
proton and the helium spectra seems to favor this model. Future 
measurements of secondary CR spectra and of the Galactic diffuse 
7-ray emission at TeV energies can provide a deeper understanding 
of the problem. 
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